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ABSTRACT 


Thermospheric airglow emissions in the wavelength range of 1800 A to 3400A and 
altitude range of 100 to 330 km were identified. These emissions compose the preliminary 
data analysis for the Naval Post Graduate School MUSTANG rocket experiment scheduled 
for launch from White Sands, New Mexico in February of 1990. Identification of 
emissions in the wavelength range of the MUSTANG experiment was accomplished using 
experimental results to support calculated emission intensities. All emission bands that will 
contribute to the experimental spectra, as well as those features that will be distinguishable, 
have been tentatively identified. These results confirm the contention that the instrument 
wavelength range allows observation of the major midultraviolet airglow emission bands. 
The extended wavelength and altitude range for this experiment allows a unique 
experimental opportunity to observe all major neutral atmospheric constituents in one 
spectrum. The primary emission band systems that will be observed are the Nz Lyman- 
Birge -Hopfield bands, the O2 Schumann-Runge bands, the NO Gamma bands and the O 
2972A line. Secondary contributions will be observed from the N2 Vegard-Kaplan and 
2nd Positive bands, the O2 Herzberg bands, and the NO Epsilon and Delta bands. These 
combined spectra will insure that the MUSTANG experimental goal of measuring N?2 and 


NO densities and determining EUV and photoelectron fluxes can be met. 


TABLE OF CONTENTS 


INTER COIBUIRC TPCQIN|eo sab coe N SOC SSN SHEO. SOG MABR (Se 0 Ree << e yt 
AS, TSC ISUS GG GS) see Z 

i TTR ORS (COTE EUS Ree enn on hos One oC etc ote tht oe 2 

J] S08 TROPIC RSI27 8 [ERC BOIS a soc onan oc oe ee ore 3 
A. THE HOMOSPHERE AND HETEROSPHERE ..................cccceeceeeeeeees 3 
Bene MlONOSPHERIC ENVIRONMENT cocccccciccccsecsccecesceesssceseeceeeeees 5 

Co AN BCE ROIN 01s)? 8 DENS RGHEN BY: ) @) Sh SSC taht eee ne 11 

il, “Tae? O31 RRRRIRGe eee eeh at. "5 Cae ene ee 14 

%, 0G ME hh coe. ree 15 

ym co MME DART CT olan, ed epee cs ke OSG Siew <nlsvs v9 See se asses ees 15 

Soy WT SPARS 2100.19 IRI la 16 

il. THE SOURCES OF THE ATMOSPHERIC AIRGLOW ................. eee Ly 
Pee DCEO SWINILIGHT SCAT PERING) | p..:5..c..ccccsccoscs.c.c0ceseccecseeesceces 19 

Sree eerrON IMPACT EXCITATION cicsecc..ccccscicccscscccsccsecsecsscescosees Pap 

Cre OCHPEMICAI “REACTIONS ................sccrercrssrscressosecsccsscnveess 22 

Pee OCW SPECTROSCOPY oi... ccciccsliceseecovencescoccngiieiacecsucsseccscseccess 24 
Pee Wr TRANSITIONS iiioici cise tds iste seis sesscncseeeececeecssosesees 24 

lee Otational Wevels sands ransitiGnS ocemeesesiisedearsisssssisnisesdsss.+.+- ZS 

Zoey iWtationalelcevels and Transitions .......<..0<02.2.ssecsisecsesecasessoee0. 28 

3. Electronic States, Nomenclature and Transitions .....................00.. 30 

SUM CONOMIG ME LALCG Bote cece nse esos Whee sees ss. ceseiecceoecsaaee 32 

DIE CILOMICe ALANGIUIOM S245) aiee seer ive rssceseesessescsdeseseens 32 

c. Electronic State Nomenclature .................cccccesccesceeeseeees 34 

al, “SEIGC INT, ROU) USER RB Bee ase oe ReRe 36 


Vv 


B. SYNTHETIC SPECTRA oie nee ee 38 


L. OVErVIEW oo... cccccccccccececcececer coucse cee ese lentetee a= 38 
2. Line PoSItioms..........20200c0cesse+++ens0 sss een erent tee 38 
3. Molecular State Populations.....-.................-.-:-:sss eee 46 


V. DAY GLOW EMISSION FEATURES IN THE 1800-3400A WAVELENGTH 


RANGE ooicceccccccscencseceeccecceueeecesensea seni esse lest ete === ===. 53 
A. MOLECULAR NITROGEN 22.2012 iscsi 54 

l. Vegard-Kaplan Bands  ........:.00.0.02...c2sesese se 37 

2. Noa Second Positive Bands ...........................1.002sseneneee 61 

3. Lyman-Birge-Hopfield Bands ..............:.2:12..---....s44.ssn 67 

B. MOLECULAR OXYGEN ......:000:1.....-::ceeeeeee ee: ee 70 

1. Schumann Runge Bands .......20sossccsssccsessosscgeieee ee eee 71 

2. Herzberg Bands. ........0is0s.6.csesss:.00-20 une ee 74 

C. NITRIC OXIDE  ....0.cc.si teense tees os seen 76 

l. Gamma Bands  ....ses.ccisosisecscsecnevaseegeeaceee eee 719 

2. Epsilon Bands. «..........csessescess+<00ctecaseces cesses eee 81 

3. Delta Bands. .....c.cccsi ee eee 83 

D. ATOMIC OXYGEN | oi.cc.cceliiesecscees- ues ece el eee eee 84 

VI. CONCLUSION | .....cciboc cece seis ccwcti cae olnteterdanenateer tee eee 87 
A. SUMMARY OF FINDINGS oon... iiccces:.s 2s ee 87 

B TOPICS FOR FURTHER INVESTIGATION ....................22 yee 88 
BIBLIOGRAPHY _....iccccccceucceciecsecccccecscees ese deee=ce eke eee eee rn 89 
INITIAL DISTRIBUTION LIST...............400.........:..0:0-00 seattle nee 93 


v1 


LIST OF TABLES 


TABLE PAGE 
2.01 IONIZING WORK FUNCTIONS FOR NEUTRAL PARTICLES.......10 
2.02 IONOSPHERIC REGIONS; PEAK HEIGHTS ........cccccceeceeee. 12 
4.01 APPROXIMATE MOLECULAR SELECTION RULES ........0.00000.. 37 
4.02 QUADRUPOLE AND MAGNETIC DIPOLE SELECTION RULES ...38 
5.01 RELATIVE POPULATIONS OF THE A2Sty STATE .....0.ccccee 58 
5.02 INTENSITY FACTORS FOR N2(C3Py - B3Pg)TRANSITION ........ 63 
5.03 O2 SR(3,7,BAND INTENSITY CALCULATIONS 0.0.0.0. 73 
5.04 INTENSITY FACTORS FOR 02(HERZBERG)TRANSITION ........ 75 
BPOSN OXYGEN (14S) SOURCES ...0.00..........cccccscssssssssssovsececccssesenes 86 
6.01 MUSTANG EXPERIMENT: OBSERVABLE EMISSIONS ............. 88 


Vil 


LIST OF FIGURES 


FIGURE PAGE 
2.01 Atmospheric Regimes .2GG0.000..ciscsceecsssasescs sess eee n= nn 4 
2.02 Generalized Ionosphere oiceci eee cscseceree eevee ieee 6 
2.03 F10.7 and Sunspot Index <..............c2232%..............:.:2-)e eee 7 
2.04 Atmospheric Temperature” Ramges........ oe 8 
2.05 Projected Launch Date Tempefature Profile..............................:.0seesse nen 2 
2.06 Solar Radiation 250-1250A...........0......2e ee 11 
2.07 Chapman Layersicsiccs.ciccc.cccsce+ contre eee eet ene ene te a 13 
2.08 Primary lonospheric Constituents-:.222-................20.0.c01es-+0.0s2 ee eee 14 
3.01 Earth Ultraviolet Airglow and Absorption/Emission Regimes......................0000+ 18 
3.02 Dissociation and Predissociation Paths...................................)3 ee 21 
4.01 Dumbbell Modeloo:...............:..2. oe sets: sercteseseee sence 25 
4.02 Vibrational-Rotational Energy Level Diagram...................cccccescesccneceeeeseeees 30 
4.03 Stable and Unstable Potential Energy Curves..................cccccescesceecceceeeeeseees 31 
4.04 Rotational, Vibrational, and Electronic Levels...................cccecececscecececeeeeees 33 
4.05 Precession about the Internuclear AXis...............:..:css:cc2se0sscerese ese see 35 
4.06 The Symmetric Top. «..........cccc:ciss.atsecueenssseesesseeneeeieeetn nen eee 42 
4.07 Hund's Case (8) oi.ccccccccscccsdcoevcssctsesee sss ccsiseeescens veel asset ete tee ae===ann 43 
4.08 Hund's Case (b)...:0iic.c. ccs. cescttedeseteee nes ccee nen eee eee 44 
4.09 A-—Type Doubling.for the 71] State....................ce.022..1- =.= 45 
4.10 Synthetic Spectrum for the (1,0) Gamma Band ......... eee e eee eeees 52 
5.01 Photocathode Efficiency..........2..0s0.0:00sesesreesesses ane dest eeeiey eet nnn 54 
5.02 Nz Potential Energy Curvee............02::i0:.000.0ssesseses esses 56 


Silo Wp Jeweress ibe Well BYE e101 iesnseseneees ss Seseeee Seen ee af 


Sena MMR Ie CaN Reap Iai VES ANLGS:--2.........0csnaceecssssavetscocsacccsssosecnseeescccasasersoscees 61 
5.05 Electron Impact Cross-Section vs. Electron Energy...................cccceeeceeeeeeeeees 64 
SOMENCCtFOMMENeEOY VS. BIECITON PIUX...............c0.c0-csescoceseccecsscnsnscscsccssesescess 65 
5.07 Electron Energy vs. (Cross-Section x Electron Intensity) 0.0.0.0... eeeeceeeeee 66 
MEME EE SCCOMG FOSITIVE ISAINGS..... cgereiaiices. +c.ssehcsersesvecseiscdeevs coedssedestvesesee. 67 
5.09 Il, - 1X, Branches for Magnetic Dipole and Quadrupole ..............:::ccceeeee 69 
Ls US 5) ye JUBTE i gn 8 V0 0 ie a a 69 
SEMME EE OLETICIANMIE MEE LY MUDIAQEANML 22.12.22, c--ccccrcesescreosesnscrseciesetssesssevescsesereses 70 
Se MING TY lee VEIN IDTAQT ANN | 2222s. ee ceca secacescceceossescceseacsocceasshvaschaccerses vvcsess 71 
PEN HSORDtLON “CTOSS-SECUION LOL OD <. 2.2... .22.ccccccecccccescs sca cecssscecsesensscscvcescees Ue 
MMMEO) Se SCHUIMANM-NUMNGE BANS... .....05--.ccccececcacccsecececdessscccrsocecsesseecsccccseres 74 
as aE Te Ty SS LENS rae ea to ne ee eae eas oe is cides widinv cw SSUES wea teeeeee 76 
SENG OtCIItial MELOY IDTARTAMY .ciie.....1.1....ssccessscssccssecsscatnescosessessescsacaess Ty 
Pepe Oemiemergy evel DiaCTaAny si e.ii 1. y.2cecccccsncocscaceceecccosscscotecccnvesesedseeesccess 78 
MMMM CO) eCereeETDITTA IS ATIGS) 5550.5 oeensceeseececese esta tn cces+-steswedcodososcetsotescsctescserceccdns 81 
TS OME TPL OTL ES ALICE errr ee) ee atte de islets + shlsa 93'ee vsicw sige ome sacnc scl dosddeleedescsesecee 82 
Seem (COME i (clings ANCL See aes cn ese acgnene es cess osacae see esbensesegedesseecsecevscens 84 
CO MMISCIET SY NEC VEN IDLAU TAN © 6. .2.ccscsccccses.-seocccstscoscdadsdsescscadsecccceserdeduceseseces 85 





I. INTRODUCTION 


There are many factors motivating the investigation of ionospheric photochemistry. 
The Department of Defense has several areas of interest in this field. A 1986 Joint Chiefs 
of Staff memorandum identified fifty critical geophysical parameters that need to be 
determined. The determination of electron density was ranked fifth on the list (MJCS, 
1986). The DoD's concern in furthering the understanding of the ionosphere is based on 
three known needs: 

- Measure the ionosphere for over the horizon radar applications. 

- Measure the ionosphere for uni-directional communications. 

- SDI beam weapons require identification of low absorption windows for 

propagation in the upper atmosphere. 
Additionally, there is a purely scientific interest in fully identifying physical processes in 
the upper atmosphere. 

These motivating factors have led to the Naval Post Graduate School rocket 
experiment called MUSTANG. The mission of this experiment is to investigate 
ionospheric photochemistry. Three specific experimental goals have been established. 
These are: 

- measure NO density in the E-region ionosphere. 

- measure N? density in the E-, Fl, F2 regions. 


- infer solar EUV flux and photoelectron flux. 


This is a joint rocket experiment containing two ultra-violet spectographs. The 
Naval Post Graduate School is responsible for a spectograph measuring ultraviolet 
radiation from 1800 angstroms (A) down to 3400A. The Naval Research Laboratory is 
responsible for a spectograph measuring from 500A down to 1500A. The launch vehicle is 
a Black-Brandt sounding rocket, which is scheduled for a February, 1990 launch date at 
White Sands, New Mexico. NASA is responsible for the launch. Observations will be 
made from 100 to 330 km in altitude. 
A. THESIS GOALS 

The primary objective for this thesis is to identify emission band systems with 
features between 1800 and 3400 A in wavelength, that are observable from 100 to 330 km 
in altitude. The accomplishment of this thesis objective wul assist in verifying that the 
MUSTANG experimental goals are attainable. 
B. THESIS OUTLINE 

Chapters II through V lead to the identification of those atmospheric constituents 
with emission features in the 1800-3400 A range. Chapter II describes the experimental 
environment which is the ionosphere between 100 and 400 km. Chapter II describes the 
general physical process that are the principal methods of transferring solar flux energy to 
create the ionosphere and airglow. Chapter IV describes two facets of molecular 
spectroscopy. The first of these is molecular transitions. The second is synthetic spectra. 
Chapter V identifies which emission band systems for the ionospheric constituents (N2, 
O72, NO, O) are observable in the experiment. Finally, Chapter VI provides conclusions 


and suggestions for further research. 


Il. THE TONOSPHERE 


The purpose of this chapter is to describe the region in space that is the focus of this 
experiment: the ionosphere. It is helpful to consider the entire atmosphere, before 
proceeding to focus on the ionosphere. The general atmospheric regimes are shown in 
Figure 2.01. Several of these regimes will be discussed in this paper. The first two 
regimes of concem are the homo- and heterospheres. The division of the atmosphere into 
these two sub-spheres is a result of significantly different processes that influence the 
layered structure of the ionosphere. Therefore, the homosphere and heterosphere will be 
discussed initially, followed by a description of the ionosphere, its sources, and its 


sublayers. 


A. THE HOMOSPHERE AND HETEROSPHERE 
The homosphere is the region of the atmosphere where turbulent mixing or eddy 


diffusion is the dominant dynamic mixing process. Atmospheric constituents N2, O2, Ar, 
He,and H, are uniformly distributed in this region. The homospheric air density falls off 


exponentially with increasing altitude as shown below: 


p (Zz) = po e-@/H), (2.1) 


The altitude is given by z and H is the atmospheric density scale height. This is defined 


below as 


(2.2) 


where kg is the Boltzman constant, T is atmospheric temperature, g is acceleration due to 
gravity, and m is the mean molecular weight of the medium. This scale height is common 
to all constituents in the homosphere. The expression defining the scale height can be 
simplified by eliminating the second term when the temperature profile becomes isothermal 


and the temperature gradient goes to zero. 
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Figure 2.01 Atmospheric Regimes (Tascione, 1988) 


The turbulent mixing processes associated with the air-earth interface diminish with 
increasing altitude. At around 100 km the dominant mixing mechanisms become molecular 
diffusion rather than eddy diffusion. As shown in Figure 2.01 the boundary is known as 
the turbopause. Atmospheric constituents begin to separate out based on molecular 
density, much as oil and water separate. The number density for each constituent falls off 
exponentially in proportion to its scale height. The density scale height is no longer based 
on mean molecular weight, but rather each individual molecular weight. This region is 
called the heterosphere. 

B. THE IONOSPHERIC ENVIRONMENT 

With increasing altitude, molecular diffusion rather than eddy diffusion becomes the 
dominant mixing process. This process combines with an increasing solar flux, and a 
decreasing atmospheric density to produce ions that have a sufficiently long lifetime to 
support the existence of a continuous plasma. This plasma is known as the ionosphere. 
The general relationship between incident solar flux and atmospheric gas concentration is 
shown in Figure 2.02. As is shown in this figure the lower limit of the ionosphere is 
typically 50 to 70 km (Tascione, 1988). The relatively wide range is primarily due to 
fluctuations in solar flux. The bold curve in Figure 2.02 represents the number density of 
ions over altitude The ionizing radiation is primarily the extreme ultra-violet and x-ray 
portion of the solar spectrum. The general shape of the ionospheric region indicates that 
ion density is negligible at very high altitudes. It then begins increasing in number density 


to a peak and decreases to a negligible intensity by an altitude of 50 km. 
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Figure 2.02 Generalized Ionosphere 

The fluctuation in the upper and lower limits of the ionosphere, as well as 
fluctuations in the ion density vary in proportion to the level of solar activity. The vast 
majority of solar flux may be treated as a black body spectrum, based on the solar surface 
temperature, of about 5700 degrees kelvin. However, the high energy end of the 
spectrum, in the ultraviolet and x-ray regions, as well as the low energy end at radio 
frequencies exhibits significant fluctuations (Torr et al, 1979 for solar flux data over the 
solar cycle). The composition and vertical limits of the ionospheric layers are significantly 
influenced by this varying solar flux. The variable portion of the solar flux varies in 
predictable cycles. This solar cycle has a period of approximately eleven years. An 
accepted indication of solar activity is the measurement of a solar radio emission whose 
wavelength is 10.7cm. This measurement is designated the F10.7 and is taken at Ottawa, 
Canada. Both the F10.7 number and the sunspot number are shown in Figure 2.03. As 
the figure indicates, the F10.7 number tracks closely with the sunspot number. 
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Figure 2.03 F10.7 and Sunspot Index(NRL, 1989) 


The variable solar flux results in atmospheric temperature fluctuations, especially in 
the upper region of the atmosphere known as the thermosphere. The thermosphere is 
primarily heated by atomic oxygen absorption of extreme ultra-violet radiation (Tascione, 
1988). Because, the extreme ultra-violet flux varies widely, the temperature of the 
thermosphere also varies. As is shown in Figure 2.04 this results in thermospheric 
maximum temperature ranges from 500 to almost 2000 degrees Kelvin. The predicted level 
of solar activity allows an atmospheric temperature profile to be estimated for a given 
location and date. An atmospheric temperature profile has been prepared for the 
MUSTANG experiment in Figure 2.05. The profile was obtained using the neutral 
atmospheric model MSIS-83 (Hedin, 1983). The model assumes a launch site latitude 33° 


N, and longitude 105° W at White Sands, New Mexico. It also uses the predicted F10.7 
value for an early February 1990 launch date of 215. The model also considers the launch 
time of day, which will be 1000 hours. The solar zenith angle is 72.47 degrees. Active 


solar conditions are used in the model. 


LOW SOLAR 
ACTIVITY 


HIGH 
SOLAR 
ACTIVITY 


ALTITUDE (km) 


THERMOSPHERE HETEROSPHERE 


MESOSPHERE 


HOMOSPHERE 


STRATOSPHERE 
=X 





TROPOSPHERE 
i) 500 1000 1500 2000 


KINETIC TEMPERATURE (°K) 


Figure 2.04 Atmospheric Temperature Ranges (Tascione 1988) 
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Figure 2.05 Projected Launch Date Temperature Profile (Hedin, 1983) 


Solar radiation of sufficient energy to ionize atmospheric constituents is the primary 
ionizing mechanism throughout the ionosphere. The far and extreme ultraviolet, portion of 
the solar spectrum contain enough energy to ionize the principal constituents of the 
ionosphere. The minimum ionizing energy requirements for the constituents are 
determined by the ionizing work function of each species. These work functions are given 


in Table 2.1 for the primary atmospheric inhabitants 


ve 


ON NG WORK N ON OR N RA tis 


Table 2.02 indicates that the ionizing radiation requirement is at the highly variable, x-ray 
and extreme ultra-violet portion of the solar spectrum. The principal ionizing portion of the 
solar spectrum is shown in Figure 2.06. It is interesting to note that the strong Lyman 
Alpha solar emission at 1215A will only meet the ionization energy requirement for the 
brelatively scarce NO molecule which results in its increased importance in especially the D 
region of the ionosphere. Photoionization is a direct energy absorption. The 


photoionization reaction is shown below for molecular oxygen: 


Oo +hv => O2T+e 
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Figure 2.06 Solar Radiation 250-1250 A (Banks & Kockarts, 1973) 


This is the initial ionization process for each of the major atmospheric constituents. The 
kinetic energy of the electrons produced may be quite significant. These electrons are 
known as, photoelectrons. They will participate in one or even several additional 
collisions, before finally reaching thermal energies of around .025 electron volts (Banks 
and Kockarts, 1973). 


C. THE IONOSPHERIC LAYERS 


A closer examination of the ionospheric zone shown in Figure 2.02 reveals that it 1s 


not a smooth curve. The vertical structure of the ionosphere is changing constantly. It also 
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varies with time of day, solar activity, as well as with latitude. Despite fluctuating 
dramatically over time, and location, the ionosphere is always present. Furthermore, the 
ionospheric plasma has four distinct regions which are primarily the result of three 
competing factors: 

- the atmospheric absorption of solar flux varies with altitude 

- recombination of ions is density dependent 

- atmospheric composition varies with altitude. 

The ionospheric layers are designated the D, E, Fl, and F2 region. The layer 
boundaries are not distinct. Although they are all present during the day, some do not 
survive throughout the night. Their general peak intensities correspond to the altitudes 
shown in Table 2.02. 

TABLE 2.02 
IONOSPHERIC REGIONS: PEAK HEIGHTS (km) 

D -90 

E -110 

Fl - 200 

F2 - 300 
These layers, are called Chapman layers after Sidney Chapman, who first modeled them in 
1931 (Chamberlain, 1987). The general shapes and locations of the Chapman layers are 
shown in Figure 2.07. The solid lines in the figure indicate electron concentrations under 


active solar conditions. The dashed lines indicate solar minimum conditions. 
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Figure 2.07 Chapman Lavers (Chamberlain, 1987) 


Before describing each ionospheric region it is necessary to describe the distribution 
of neutral species in the heterosphere, which overlays the E and both F layers of the 
ionosphere. The relative abundance of the neutral particles is distributed differently in each 
layer. The neutral atmospheric particle distribution is shown in Figure 2.08. Below about 
175 km N?2 1s the dominant species present. Above this altitude atomic Oxygen becomes the 
dominant species. Oz and NO remain relatively minor constituents although as will be 


shown in Chapter V, both have significant roles to play in ionospheric photochemical 


processes. 
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Figure 2.08 Primary Ionospheric Constituents 


(Hedin, 1983 and Cleary. 1986) 


The dominant processes that produce each 1onospheric region is described below. 
A more detailed treatment of each region may be found in [he Theory of Planetary 
Atmospheres, by Chamberlain (1987). 

1. The D Region 

The D region extends from 50 to about 90 km, with its peak near 90 km. 

The relatively high number density in this region results in a short mean free path between 
neutral particles, and electrons. The daytime electron density is about 1.4x104 electrons 
per cubic centimeter (el cm-3). This produces a high recombination coefficient that causes 


the region to disappear at night. Solar Lyman-alpha radiation deposited in the region ts the 
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principal ionization source for the NO molecule. Although photoionization is the dominant 
process, other interactions contribute to ion density. These include primarily three-body 


attachment as shown below: 


Or +e+M> 07°+M. 


A reaction of secondary significance is radiative attachment such as the example shown 


below: 


Or+e => Or + hv. 


Principal loss mechanisms include photodetachment, which is the inverse of radiative 
attachment; collisional detachment, the inverse of three-body attachment; and mutual 
neutralization, such as Xt + Y- > X + Y. This last process indicates that the D region is 
characterized by a significant number of negative ions as well as electrons. 
iss The E region 
The E region extends from 90 to about 130 km. The nominal peak height is 
110 km. The daytime density of 1.5 x10° (el cm-3) decays down to a night time density of 
1x104 (el cm-3) or less. The decreasing recombination coefficient allows its nighttime 
survival. Ionization primarily is the result of photoionization by solar extreme ultraviolet 
and soft x-rays. The ionization of atomic oxygen becomes significant. Negative ions are 
no longer significant contributors in the region. The principal ion loss mechanism is 
dissociative recombination such as: 
NOt +e=N+0O 
o.. Fl region 
The Fl region extends from 140 to 200 km. Its daytime electron density is 


greater than the E region (2.5x10° el cm-3vs. 1.5x10> el cm-3) However, the principal 
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ions formed are OF and NOt which require an ionizing solar flux from the Lyman 
continuum as well as He emission lines. These sources disappear with sunset. Due to the 
relatively short ion lifetimes, the F1 layer disappears quickly at night. O* readily charge 
exchanges with O2 and NO. These exchanges are given below: 
Ot + Op > O+ O2t 
Ot + No => NOt +N. 
Once in these molecular forms the principal loss method is once again dissociative 
recombination. 
4. F2 region 

Beginning at 160 km, the F2 region extends all the way up and beyond 600 
km, with a nominal peak density at about 300 km. Its day time density is around 109 (el 
cm-3). The night time density drops by about one order of magnitude, to about 109 (el cm- 
3). In the F2 region three significant factors adjust the ionization profile. First, the 
atmosphere becomes optically thin to most ionizing radiation; which adjusts the ion 
production equation. Second, 1on-atom interchange becomes the dominant loss process 
instead of dissociative recombination. This results in the loss term diminishing as 
atmospheric density decreases. Ion production is from the same processes as the F1 layer. 
Atomic oxygen is the principal constituent and ion source in the region. However, since 
the loss term is diminishing with altitude the ion number density would increase indefinitely 
with height except for the effect of the third factor; ambipolar diffusion. The F2 region 
reaches the point where the chemical lifetime becomes longer than the diffusion lifetime. 
As a result positive ions diffuse downwards out of the region. This limits ion density 


towards the upper end of the region. 
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Il. THE SOURCES OF THE ATMOSPHERIC AIRGLOW 


x BACKGROUND 

The MUSTANG experimental objective to investigate ionospheric photochemistry 
is best accomplished by observing radiative emissions in the earth's airglow. Airglow 
emissions are omni-directional. An instrument's viewing angle is oriented to observe 
emission features without being overwhelmed by the solar flux. Specific emission features 
can be observed to identify the atomic or molecular transitions and therefore the molecular 
density of the species. Airglow emissions can then be converted to density profiles. These 
number densities can then be used to infer the electron and ion densities. 

The sum of all radiative emissions from the extreme ultraviolet to the infrared 
comprises the atmospheric airglow. The airglow is optical radiation that arises from several 
types of interactions. Airglow emissions have three general sources: 

- direct sunlight scattering 
- electron impact excitation 
- photochemistry of neutral constituents. 
The primary source is direct sunlight scattering, or photoexcitation. Electron impact 
excitation is of secondary importance for allowed transitions, but it is the primary excitation 
mechanism for forbidden transitions. Some radiation is the result of a combination of these 
factors. 

Solar ultraviolet radiation is primarily responsible for the airglow to be observed in 
the MUSTANG experiment. Figure 3.01 shows a typical nadir spectrum of the earth's 
airglow from 500 A to 4000 A. The primary EUV and FUV emission features are shown in 
this figure.The strong O2 and O3 absorptive bands shown in the figure result in the 
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LOG INTENSITY (Rayleigh/10 A) 


suppressed airglow at wavelengths corresponding to these absorption bands. The figure 
also illustrates the emission bands to be observed by the Naval Research Laboratory 
instrument . The strong Lyman Alpha, and atomic Oxygen emission features are clearly 
visible between 500 and 1500A. Conversely the NO and N> features in in the 2000 to 
3000 A range are not visible in a nadir view. This is an illustration of why the experiment 


must look across the limb to observe emissions in the 1800 to 3400 A region. 
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Figure 3.01 Earth Ultraviolet Airglow and Absorption/Emission Regions 


The airglow is distinguished from another emission phenomenon, the aurora, by 
the energy source of the aurora. The energy sources for the aurora are solar electrons and 
protons, as opposed to solar flux as the primary source of the airglow. This results in the 
auroras confinement to polar regions and its sporadic occurrence. Conversely the airglow 


occurs continuously, over the entire sunlit portion of the atmosphere. 


A. DIRECT SUNLIGHT SCATTERING 

The primary process that results in the airglow is direct sunlight scattering, or 
photoexcitation. This process is similar to photoionization, the major ionization source. 
Photoexcitation occurs when the incoming photon energy matches the quantum energy 
levels required for a given atom or molecule. The photon will then be absorbed and leave 
the particle in an excited state. The normal lifetime of a particle in an excited state is very 
short; about 10-8 seconds (Eisberg, 1985). 

The particle in the excited state can return to its desired ground state by emitting one 
or several photons. Photoexcitation that leads to the emission of photons is known as 
fluorescence. The entire process is known as fluorescent scattering. A special case of 
fluorescent scattering is resonant scattering. In the case of resonant scattering, the 
excitation and emission paths are the inverse of each other. A single photon, of the same 
wavelength as the excitation photon, is reradiated. The particle retum path in general, may 
be along a direct or an indirect path. Multiple photons are emitted, for the general case 
where intermediate states lie between the excited state and ground. Only electronic 
excitation levels are available to atoms. Molecular excitation may result in electronic, 
vibrational, or rotational transitions, or combinations of these three. Energy requirements 


vary for each type of transition. Electronic transitions require between one and ten electron 
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volts, while vibrational transitions require between one tenth and one electron volt, and 
rotational transitions are a factor of ten less than vibrational transitions. 

Transitions are considered allowed or forbidden depending on their compliance 
with quantum selection rules. The actual selection rules that apply to molecular transitions 
will be discussed further in Chapter IV. A transition that meets all selection rule 
requirements is an allowed transition. A transition that violates a selection rule at one level 
may be forbidden; but may or may not be strictly forbidden. For example a transition 
forbidden as electric dipole radiation may be allowed with a lower probability as a magnetic 
dipole, or electric quadrupole transition. 

An excited state whose only path to the ground state is through a forbidden 
transition is called a metastable state. Usually, the particle can still reach the lower state, 
but the probability of doing so becomes much smaller than if the lower state was an 
allowed state. Lifetimes in the excited state can become relatively long. They can be on the 
order of seconds as 1s the case for one of the emission bands to be examined in Chapter V. 
Metastable states are much more likely to transfer energy as a result of a collision with 
another particle prior to emission. This process is called quenching. 

It is also of note that molecules in an excited state may dissociate, or predissociate, 
rather than emitting a photon. An absorption leading to an excited state above dissociation 
or predissociation energy, will dissociate rather than emit a photon. These processes are 
best shown graphically. Figure 3.02 shows two simplified potential energy diagrams for 
the Oz and NO molecules. As shown in the figure, an O2 molecule excited to the A3Z+, 
state at a vibrational level above the dissociation threshold energy will dissociate to the 
OP) + OP) atoms. This may be viewed as a path to climb out of the electronic potential 
well. Conversely, a molecule of NO excited to a vibrational level above the 


predissociation energy of the C2IT state may predissociate along the a“IT path to become 
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N(4S°) + OP) atoms. This may be viewed as a hole in the potential energy well. The 


result of both of these processes is that absorptions will not result in emissions. 


O2 Potential Energy Diagram 
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Figure 3.02 Dissociation and Predissociation Paths (Chamberlain, 1987) 
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B. ELECTRON IMPACT EXCITATION 

The primary excitation source for atoms or molecules is solar radiation. However, 
photoexcitation is forbidden for some of the constituents in the upper atmosphere. In this 
case, electron impact can cause the excitation of the particle to an excited state where it can 
than emit a photon upon undergoing a forbidden transition. Electron impact scattering is an 
inelastic collisional process. The probability of an emission resulting from an electron 
collision with a neutral particle depends on several factors. Electron impact excitation ts 
dependent on the molecular number density, electron impact cross-section, and the electron 
flux. Since the electron flux is not monoenergetic, and since the molecular electron 
absorption cross-section is dependent on electron energy, an actual calculation requires an 
integration over available energies. This will be done for several forbidden emissions in 


Chapter V. 


Cs PHOTOCHEMICAL REACTIONS 

A third relatively minor source of airglow emissions is radiation resulting from the 
photochemistry of neutral constituents. Photochemical processes are initiated in the same 
manner as photoabsorptive processes. A photon is absorbed leaving the absorbing 
molecule in an excited state. This excited state could be an excited rotational or vibrational 
level or an excited electronic state. It could also be a dissociated or ionized state. In this 
excited state the molecule may enter into secondary chemical reactions that would not be 
available from the ground state. Upon completion of the chemical interaction the new 
species emits a photon. This process is called chemiluminescence. 

Photochemical processes are complicated because they are non-local 


thermodynamic equilibrium processes. Each rate of reaction must be calculated separately. 


ZZ 


Addressing each photochemical process is beyond the scope of this thesis. 
Chemiluminescence is an insignificant contributor to the day airglow. It becomes important 


at night when other excitation sources are not present. 
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IV MOLECULAR SECTROSCOPY 


Chapter IV addresses molecular transitions, and synthetic spectra. These are two 
broad regions that comprise molecular spectroscopy. The vast majority of emission 
processes to be observed in the MUSTANG experiment are a result of molecular 
transitions. Furthermore, this chapter is exclusively concerned with diatomic molecules. 
Diatomic molecules are the dominant type found in the atmosphere. Molecules other than 
diatomic introduce coupling schemes beyond the scope of this thesis. More specifically, 
emissions will primarily occur as a result of vibrational transitions between molecular 
electronic states. This type of emission process 1s one of several possible types of diatomic 
molecular transitions. Molecular transitions occur between rotational and vibrational levels 
as well as between electronic molecular states. It may also be between some combination 
of these three. Each type of transition is characterized by quantum energy jumps rather 
than a continuous emission spectrum. All three types of molecular transitions will be 
addressed, with emphasis on vibrational transitions between electronic states. Additionally 
the compilation of synthetic spectra will be addressed. A more detailed account of 
molecular transitions is provided by Herzberg, in his definitive volume, Spectra of 
Diatomic Molecules (Herzberg, 1950). 


Aas MOLECULAR TRANSITIONS 

This section is concerned with rotational and vibrational levels, as well as electronic 
states. It is also concerned with transitions between these states. These transitions range in 
energy from the least energetic rotational transition, to the most energetic electronic 


transition. They are in effect a nested sequence of transitions. A pure rotational transition 
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can occur between two rotational levels in the same vibrational level and electronic state. 
However, as vibrational and electronic transitions come into the picture, subordinate level 
transitions must be taken into account. 
1. Rotational Levels and Transitions 

A diatomic molecule as shown in Figure 4.01, can be represented as a 
dumbbell. The two masses m, and m2 are connected by a rigid rod of length r with the 
center of mass at s. The two masses represent the atomic nuclei, where the vast bulk of the 
atomic mass is concentrated. As opposed to an atom, the molecule can undergo two forms 
of motion. It can rotate around the center of mass axis and the two nuclei can vibrate 
relative to one another along the internuclear axis. The two masses are not in fact ngidly 
bound. This perturbation to the rigid rotator and vibrator will be addressed later in this 
chapter. At this point it is adequate to note that the internuclear axis will vary with the 


speed of rotation, or strength of oscillation. 





Dumbbell Model of a Diatomic Molecule 


Figure 4.01 Dumbbell Model (Herzberg, 1950) 


The simplest form of motion the two components can undergo is rotation 


around the center of mass axis S. Quantum theory requires the velocity of this rotation to 


ZS 


take on discrete values. Isolating the molecular motion to pure rotation results in the 


system's angular momentum as 
P=Io. (4.1) 


The classical kinetic energy of rotation would be E = 5 Iw. The classical energy becomes 


the expression 


P is from equation 4.1, and the moment of inertia I, is dependent on the reduced mass. 


The solution to the Schrodinger equation yields available energies of the form shown 


2 
p+ J(J + 1) 
8nx2 I 


(4,3) 

In Equation 4.3, h is Planck's constant, I remains the moment of inertia, and the rotational 
quantum number J is introduced. J is limited to positive integer values. The energy of the 
system E, will increase quadratically with increasing J as is indicated in Equation 4.3. The 


angular momemtum of the system in some quantum state J becomes 


p= , NIG See a 7. Ga 


Te 


Classical electrodynamics requires a changing dipole moment to produce an emission. Any 


diatomic molecule composed of two heteronuclear atoms will possess a permanent dipole. 
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Diatomic molecules with a permanent dipole moment will produce an emission when 
rotated. The quantum mechanical system indicates that the system will have a series of 
discrete energy levels whose energy increases quadratically as J increases. A transition 
between two of these available energy states will result in a discrete rotational emission 
feature. The wave number v, is introduced at this time. Wave number is equal to the 
reciprocal of the wavelength (1/A) or frequency divided by the speed of light. Wave 
number carries units of cm-! (number of waves per centimeter). The wave number for the 


rigid rotator is given by 


The upper and lower energy states are indicated by E' and E", respectively, while h is 
Plancks constant and c is the speed of light. The rotational term F(J), is also introduced, 
and depends on the physical characteristics of the molecule and dependent on J. The 
relationship is 


F(J) = - -~ BJ +1). (4.6) 


The constant B is known as the rotational constant of each molecule and is given by 


~ Sx2el— 





Only certain rotational transitions are allowed for a given excited state J’. 
The available states are governed by what are known as selection rules. Selection rules for 


the rotational quantum number J, are AJ = + 1. This requires that J vary by one higher or 


one lower quantum level in either absorption or emission for a rotational transition. 
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Pa. Vibrational Levels and transitions 

A perturbed, diatomic molecule can also undergo vibratory motion. The 
configuration for the molecule shown in Figure 4.01 is still appropriate, however a second 
mode of motion is now considered. A classical picture would consider the two molecules 
as connected by a spring. Perturbing the system would result in simple harmonic motion. 
The equation of motion is a straight forward result of Hooke's Law: F = -kx, where k is 


the force constant . The vibrational frequency of the molecule is defined as Vogc = a = 
x as 


where | is the reduced mass. For a given vibrational state, the vibrational energy level, 


and the resultant vibrational term G(v) is given by 


Gw) = BM owed. (47) 





The factor @, equal to Vosc¢/c, is introduced in equation 4.7, and 1s a constant for each 
molecule. The vibrational quantum number v is also introduced. The vibrational quantum 
number is restricted to positive integers. It is of note that when v is equal to zero; the 
lowest vibrational state, the vibrational energy equals 5 @, rather than zero. Additionally 
the linear dependence on v, results in equidistant energy levels, for increasing vibrational 
levels. 

Vibrational selection rules, restrict available transitions between higher and 
lower states to unity; Av = + 1. Homonuclear molecules such as N2, have a zero dipole 
moment. Vibrational transitions can only occur as quadrupole transitions for these 
molecules. Quadrupole radiation has a very small probability with respect to dipole 
radiation, and is negligible under normal circumstances. The selection rule Av = + 1 is 


not ngidly adhered to for asymmetric molecules such as NO. Dipole allowed transitions 
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have a weak probability of emission when Av > + 1. This is as opposed to the rotational 
case where the selection rule is rigidly adhered to. 

The molecule can undergo simultaneous vibrational, and rotational motion. 
The total energy in the rotating oscillator case is approximately T= G(v) + F(J). The period 
of oscillation is on the order of 10-!4 seconds, while the period of rotation is about 10-!2 
seconds. This means that the rotational position changes slowly compared to the oscillating 
motion. For a given vibrational energy level, there are a series of rotational levels 
available. The selection rule remains AJ = + 1 for a specific vibrational level. The series 
of vibration-rotation transitions is called a band. The AJ = 1 case forms one branch, 
known as the R-branch and the AJ = -1 case forms the P-branch. An illustrative example 


of a vibrational-rotational band system, with P and R branches is shown in Figure 4.02. 
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Figure 4.02 Vibrational-Rotational Energy Level Diagram 
(Herzberg, 1950) 
oie Electronic States, Nomenclature, and Transitions 
a. Electronic States 

The nuclei of a molecule repel each other. The molecule is held 
together by what is in effect an entire molecular group of electrons. The electron cloud 
interaction with the nuclei results in a strong axial electric field. The orbital configuration 
of this electron cloud determines the electronic energy of the system. The electrons are 
moving much faster than the nuclei, so the electronic energy corresponds to each temporary 
position of the nuclei. Work must be done against the Coulomb repulsion of the nuclei, 
and work is required to change the electronic energy, in order for the position of the nuclear 
elements to change. There must be some minimum value in order for the electronic state to 


be stable. The shape of the potential energy curves for electronic states indicates whether 
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the state is stable or unstable. Without a minimum in the potential energy curve, the two 
atoms will repel each other at any internuclear distance. This is an unstable state. Figure 
4.03 ulustrates several potential curves for different electron states for a Liz molecule. As 
is shown in the figure, the 'T]y ,!£+y and !Z*, states all possess minimum points, and are 


stable states. The 3£+, state has no minimum and is an unstable state. 
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| Figure 4.03 Stable and Unstable Potential Energy Curves 
(Herzberg, 1950) 
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b. Electronic Transitions 
In a similar manner in which purely rotational transitions were 
extended to include vibrational-rotational transitions, vibrational-rotational transitions can 
now be extended to include transitions between electronic states. The total energy of the 
molecule in a given energy state is defined as the sum of the rotational, vibrational and 


electronic components. This definition is given by 


E=Ee+E,+E,. (4.8) 


The change in energy from one specific state to another is the sum of the energy changes of 


each component as 


AE = AE. + (E,' - Ey") + (E,'- E,"). (4.9) 


For any given change in electronic state, there is a wide range of possible vibrational and 
rotational transitions. Thus, the transitions become a band rather than a single line. Figure 
4.04 illustrates the concept of nested rotational, vibrational and electronic levels. This 
figure illustrates the nesting concept between electronic states. For a given transition 
between two electronic states, there are numerous vibrational transitions possible. Each 
vibrational transition is in fact a band of all possible rotational lines for each particular 
upper and lower vibrational level. The sum of these bands form the band system 


representing emission or absorption features between two electronic states. 
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Figure 4.04 Rotational, Vibrational, and Electronic Levels 


(Eisberg, 1985) 
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Cs Electronic State Nomenclature 

The intemal molecular electric field is symmetric about the 
intemuclear axis. The molecular electron's motion in this field results in a total orbital 
angular momentum, L. As shown in Figure 4.05, L is not a constant of the motion as it 
precesses around the axis of symmetry. However, the component of L along the 
internuclear axis is a constant of the nutation. This component designated A is a better 
indication of the molecular electronic state than L. The electronic configurations, and 
indeed the labeling of successive electronic states, is analogous to the atomic case. The 
symbol A =! My]! = 0,1,2,... L is used to designate a particular electronic state. 
Furthermore, each electronic state is given a Greek letter to distinguish the molecular 


electronic state from its atomic equivalent. The lower four molecular states are: 


0-2 
1 - Tl 
2-A 
3-@. 


It should be noted that for My > 0, electronic states will exhibit double degeneracy since for 


a given A, My can be positive, or negative. 
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Orbital Angular Momentum Precession 
about the internuclear axis r 


Figure 4.05 Precession about the Internuclear axis (Herzberg, 1950) 


Electron spin in molecules results in a multiplet structure for either a 
particular band or an individual line. Just as in the atomic case, an odd or even number of 
total electrons results in a half integer or integer total spin, respectively. The resultant spin 
S is unaffected by the electrostatic field. It will precess in the presence of a magnetic field. 
When A is greater than zero, the orbital motion of the electrons results in the generation of 
an internal magnetic field in the direction of the internuclear axis. The spin interaction, S, 
with this magnetic field will precess in much the same manner as L. Also in much the same 
manner as L, a constant component Msg will exist for a given spin precession. Mg has 
specific values designated 2, where 

re Eee ae eee 
This is a total of 25+1 possible values. This also indicates the multiplicity of the state. Ina 
manner similar to atomic classification, the spin is written as a left superscript. Where A # 
0, the sum A + & is written as a right subscript. 
Some additional points of classification for electronic states are 


noted at this tme. The ground state of a diatomic molecule is labeled with X, while 
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successive excited states, of the same multiplicity, are labeled A, B, C, and so on. 
Similarly, successive excited states where multiplicity differs by two are designated by a, 
Deere. 

The nomenclature being discussed above, is repeatedly used in the 
potential energy diagrams in Chapter V where various emission band systems are 
investigated. Emission band systems are generally named for their discoverer. For 
example the N2 Vegard-Kaplan bands were discovered by Vegard and Kaplan (JPCRD, 
1977). These bands are designated as eS ME. This indicates a transition from 
the first excited electronic state A, to the ground state X for the case where both upper and 
lower states are in the 2 electronic configuration. Since A = 0 for the & state, the right 
subscript of A + 2 provides no additional information. In this case the right subscript 
becomes u and g which indicates an odd to even electronic state.The spin goes from 2 
equal to three to 2 equal to one. This violates one of the selection rules that govern 
electronic transitions. These selection rules are discussed in the next section. 

d. Selection Rules 

Selection rules governing electronic transitions are not as simple as 
those for rotational or vibrational levels. The selection rules that govern allowed transitions 
are shown in Table 4.01. Allowed transitions will comply with these selection rules. This 
will result in electric dipole radiation. The AJ selection rule gives rise to three branches 
although more branches may result from a particular transition. Depending on the 
multiplicity of the upper and lower states the three branches are designated R, P and Q. 
The + — selection rule in Table 4.01 refers to positive and negative terms. Only 
positive to negative transitions and vice versa are allowed. The s & s and a > a selection 
rules refer to symmetric and antisymmetric terms. Only symmetric to symmetric or 


antisymmetric to antisymmetric transitions are allowed. The g & u selection rule requires 
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only even states to odd and vice versa are allowed. These selection rules are subject to the 
limitations arising from compliance with coupling schemes to be covered in the next 


section (Hund's case (a) and (b)). 


TABLE 4,01 
-APPROXIMATE MOLECULAR SELECTION RULES 
(Herzberg, 1950) 
AS =0 
+o- 
Ss 
aoa 
gu 
Ae oe 
AJ=0,+1 


There are exceptions to these selection rules. One exception of 
significance is that £* to LX is not allowed but &* to 2+, and X- to 2- is allowed. Other 
exceptions will be dealt with on an individual basis in Chapter V. 

Violation of one of the selection rules in Table 4.01 does not mean 
the transition is rigidly forbidden. The transition may be allowed as magnetic dipole or 
electric quadrupole radiation which introduces additional selection rules. These mules are 
supplemental exceptions to the selection rules of Table 4.01. The Herzberg bands for 
molecular oxygen as well as the Vegard-Kaplan and Lyman-Birge-Hopfield bands for 


molecular nitrogen all violate at least one of the selection rules listed in Table 4.01 (see 
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Chapter V). Table 4.02 lists selection rules that apply for magnetic dipole and electric 
quadrupole radiation. 


TABLE 4.02 
‘VMIENTAL SELECTION RULES FOR ELECTRIC ¢ RUE 
AND MAGNETIC DIPOLE RADIATION 
(Herzberg, 1950) 
+> + 
=e 
E— sg 


ua u 


B. SYNTHETIC SPECTRA 
a Overview 
The previous section defined molecular transitions along with the 
nomenclature and selection rules that apply to molecular transitions. In this section the 
wavelengths and relative intensities of these transitions will be identified. The synthetic 
spectra that are produced by these techniques can then be used to link observed emissions 
with the emitting molecule. Unique emission features allow observers to identify the 
"fingerprints' of each emitting species. The techniques described in this section can be 
used by computer modelers to develop synthetic spectra for each emission band that will be 
observed. 
Ze Line Positions 
The magnitude of the energy difference shown in equation 4.9, divided by 


hc, for either an emission or absorption, is known as the line position. Line position 
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carries units of cm-!. Another descriptive name is the wave number. The line position o, 
follows from the change in the appropriate parameters for the electronic, vibrational, and 
rotational contributions to the total. These are: the total electronic energy Te, the vibrational 
energy F(J), and the rotational energy G(V). The vibrational and rotational energies remain 
as were discussed earlier in Chapter IV. For a given transition from an upper state 


(Te',v' J) to a lower state (Te",v' J") , the line position o is expressed as is 


o=T'- T" + Gv) - G"(v) - F'(J) - F'(J). (4.10) 


The electronic portion of the entire band system Og = T' - T" remains constant. This can be 
used to simplify Equation 4.10. The line position for a particular band system can be 


written as 


6 =0e + G'(v') - G"(v") - F'(J') - F"(J"). (4.11) 


For any v to v’ transition, there may be hundreds of rotational transitions. The wave 


number for any one of these lines becomes as is shown by:. 


Oj'j" = Oot F'(J') . F"(J"). (4.12) 


The 69 term is known as the band origin or band head. Its value is given by 


Oo = Get G'(v') + G"(v"). (4.13) 


39 


At this point the rigid rotator and simple harmonic oscillator approximations 
need to be addressed. The constants used to determine line position are given by Huber 
and Herzberg (1979). However, these constants include anharmonicity and non-rigid 


rotator effects. In this case the vibrational term becomes 


G(v) = W e( V+ ,) : WeXe(V + ), (4.14) 


The values for We and Xe are constants of each molecule and the second term accounts for 
anharmonicity. The rotational term becomes 


F(J) = ByJ (JJ + 1) - Dy J20 +1). (4.15) 


where the second term accounts for the non-rigid rotator. The factors By and D, depend on 


the vibrational level. The expression defining By is given by 


By = Be - Ge( Vv + 5) + higher order terms. (4.16) 


Similarly the expression for Dy is 


Dy = De- Be( v + >) + higher order terms. (4.17) 


The constants Be, &, De, and B are all found in Huber and Herzberg. 
The spin 2 # 0 and orbital angular momentum A # 0 contribute to molecular 


degeneracy through several coupling schemes. The rotational term values in these cases are 
highly dependent on the coupling between A and Z. In 1927, Hund (1927) first 


considered five possible cases of coupling between A and 2. The first two of these are 
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applicable to molecules with emission bands between 1800 and 3400 A and will be 
examined further. The total electronic angular momentum vector Q is the algebraic sum of 


the A and £ components. The algebraic rather than vector sum shown in equation 4.18 is 


appropriate since both terms are aligned along the internuclear axis. 


Q=/L+A] (4.18) 


Hund's Case (a) is appropriate when strong coupling exists between the 
electronic orbital and spin motion, and the intemuclear axis. At the same time electron 
motion is weakly coupled to the rotation of the molecule. This situation is completely 
analogous to a symmetric top. The electron motion around the two nuclei results in a non- 
zero moment of inertia about the internuclear axis. The molecule will precess about the 
total angular momentum vector J, just as the symmetric top does. The symmetric top is 
shown in Figure 4.06. The curved arrow in the figure represents the rotation of the entire 
molecular system about the J vector. A indicates the orbital angular momentum of the 
electrons about the internuclear axis. N is the component of the total angular momentum 


perpendicular to the internuclear axis. 
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Vector Diagram for a symetric top 


Figure 4.06 The Symmetric Top (Herzberg, 1950) 


The system in Hund's Case (a) has been complicated by the interaction of 
the spin and orbital electronic motion with the nuclear rotation and with the internuclear 
line. The !X ground state, where spin and orbital angular momentum go to zero, will be 
of the same configuration as the symmetric top. Figure 4.07 illustrates the coupling in 
Hund's case (a). In case (a) the nuclear rotation and electronic motion coupling is weak, 
while the coupling between electronic motion and the internuclear line is strong. The 
magnitude and direction of the J vector is constant, while the electronic angular momentum 
(2 and the angular momentum of nuclear rotation, N, couple to produce the nutation about 
the total angular momentum J. Additionally L and S are precessing much faster than the 
precession about J. The bold arrow in Figure 4.07 indicates the nutation of the system 
about J while the ovals at the end of the A and Z vectors indicate the more rapid precession 


of L and S about the internuclear line. 


42 





Hund's Case (a) 


Figure 4.07 Hund's Case (a) (Herzberg, 1950) 
Hund's case (b) has a substantially different coupling scheme. The spin S 
is only weakly coupled to the internuclear axis. The orbital part of the electronic motion is 
still strongly coupled to the axis but more moderately coupled to the rotational motion. 


Additionally, when A => 0 case (a) can no longer apply. The K vector encompassing the 


total angular momentum apart from spin is then introduced. K is a quantum number with 


integral values as defined below: 


K=A, A+i,A+2,.. (4.19) 


The total angular momentum J now follows from 4.19 as 


J=K+S,K+S-1,K+S-2,...|]K- S|. (4.20) 
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The system configuration and relative coupling for Hund's Case (b) are given in Figure 


4.08. 





Hund's Case (b) 


Figure 4.08 Hund's Case (b) (Herzberg, 1950) 


The amount of splitting of the energy levels of a multiplet state depends on 
the molecule's adherence to Hund's case (a) or case (b). A general rule is that multiplet 
splitting is large for Hund's case (a) and approaches zero for Hund's case (b). The amount 


of splitting can be accounted for by the expression 


Ti = Tg + AAD eae 


In equation 4.21 A and 2 are as have already been defined and A is a molecular constant 
for the particular muliplet term. 

The Hund's coupling cases discussed above are idealized limiting cases. In 
general coupling does not neatly comply with one case or another. As rotational energy 
levels increase, a transition between case (a) and (b) is observed. This is known as 


uncoupling and there are two types of uncoupling mechanisms between case (a) and (b). 


These are 
A-type doubling and spin-type uncoupling. Both mechanisms will cause an energy level 
splitting. A-type doubling occurs for Hund's case (a) and (b), for A # 0 states (II, A, ® 


etc.). This doubling is very small relative to the other energy differences involved in a 
transition. Figure 4.09 illustrates A-type doubling for two I states. The splitting steadily 


increases as J increases. 


i 


0 5 ome 7 
Figure 4.09 A-type Doubling for the 2II State (Herzberg, 1950) 
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Spin-type doubling depends on the spin S for the molecule. As rotational speeds 
increase, the rotational speed becomes comparable to the precession of S about L. Hund's 
Case (a) is dominant for small rotations and splitting is large. For large rotation the 
situation approaches Hund's case (b) and splitting is reduced. The result of this larger 
splitting due to spin-type uncoupling results in different rotational term values for different 
spin sub-states. 

a Molecular State Populations 

The calculation of the population of a specific molecular state requires 
accounting for all factors that contribute to the population. These factors vary between 
Hund's case (a) and (b). A state in Hund's case (a) is defined by n, v, and 2, while case 
(b) is defined by n, v, and N. For those states where A>0, A-type doubling introduces a 
parity quantum number p that must be accounted for. Although line positions are not 
sensitive to A-type doubling, the calculation of the the population rate and transition 
probability requires the two sub-states be considered. This is because each A-doubled 
component is individually populated. Transitions from each upper A-doubled state 
constitute a fully resolved line. 

The volume emission rate can be calculated by considering the fractional 
population of the upper state multiplied by a spontaneous emission coefficient between the 


two states. This is shown for Hund's case (b) in the equation below: 


E(n'v'N'J'p'n'v"N"' J" p") - Ni(n'v'N'J'p') A(n'v'N'J'p',n'v"N"J"p"). 
(4522) 


where E is the volume emission, N is the fractional population, and A is the Einstein 


coefficient for spontaneous emission between two states. The factors enclosed are 
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parameters defining the electronic state n, the vibrational level v, the rotational level N, the 


vector J, and the parity p. For Hund's case (a) the parameter 2 would replace the rotational 
level N. 
The number of molecules in a specific n' and v' level may be accounted for 


by the production rate divided by the loss rate in the steady state using 


NG ge NS aga) 


aviv" 


Now Equation 4.22 can be multiplied by | in the form of N(n'v')/N(n'v’) and expressed as 


TLItnst rtaat tt te Agih ree thy tavit 'v'N'J'p' yy" 
E(n'v'N'J'pin"v"N'J"p") = P(n'v' eo es . (4,24) 
dAyty" 
Vv 


In equation 4.24 the factor Ayy” is equal to A(n'v'N'J'p',n"v"N" J" p"). 

The production mechanism leading to emission for allowed transitions 1s 
photoexcitation. The production rate P(n'v'), is proportional to the solar flux and 
effective cross section for the transition to a specific electronic state. It is assumed that the 
upper state rotational population is statistically similar to that of the ground state. The 
branching ratio into various vibrational levels are given by the Frank-Condon factors. 
These factors arise from the Frank-Condon principle. A brief description of the principle 
follows. The line strength for a particular transition is a function of the relative positions of 
the wave functions between the upper and lower states. The electronic transition time for a 
molecule is on the order of 10-!© seconds. Because the vibrational transition time of 10-!4 
seconds is relatively long the internuclear distance will remain virtually constant during the 


electronic transition. Referring back to Figure 4.04, this transition may be expressed as a 
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vertical line between electronic states. The maximum probability density in the excited state 
will determine the upper state location. A strong wave function spatial overlap will 
determine the most likely lower eigenstate. The most favored location is generally at the 
end of an oscillation where the eigenfunction is a maximum. The relative strength of a 
transition from a vibrational level (v') of an excited electronic state to the vibrational level 
(v") of a lower electronic state is then based on the locations of the upper and lower levels, 
and the probability amplitude at the respective upper and lower locations. This concept, 1s 
known as the Frank-Condon principle (Eisberg, 1985). 


The production rate for an allowed transition can be written as 


aes Tez P 
Pin waa K Fog ov' fxn' Qov' (4.25) 


where fF is the solar flux (photons cm-2 sec"!A-!), 2 is the wavelength of incoming 
radiation, fxn’ is the oscillator strength for the electronic transition from the ground state to 
the n' state, and doy' is the Frank-Condon factor between the v=0 level of the ground state 
and the v' level of the upper n' state. This production equation makes the valid assumption 


that virtually all molecules are in the ground state prior to absorbing incident radiation. 


Tatum (1967) develops a detailed description to determine the relative 
populations of molecular energy levels. The development is summarized here. The 
fractional populatio art = é of a given state generally obeys a Boltzman 


distribution. The Boltzman distribution for a state can be expressed as 


: . e-E./KT 
Nis = BSE 4.26) 
YS gje"Es/KT 
j 
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In Equation 4.26, gj is the statistical weight of the ith group and the denominator is the 
associated partition function. Tatum applies this distribution to the different molecular sub- 
states corresponding to appropriate Hund's coupling cases. As an example for Hund's 
case (b) in states where A # 0, the fractional population may be expressed as 


N(niviN'J'p') _ he By g (J + 1) hceF(N) 
Nin'v') = kr PQS a1) &xPC4p)- (4.27) 


The factor ® introduced in Equation 4.27 is a constant for hetero- or homo- nuclear 
molecules. The various other factors in equation 4.27 have all been defined previously. A 


similar expression is obtained for Hund's case (b). 


The remaining factor in equation 4.24 ee) eal , 1s yet to be determined. 


Save 
The factor Ayy" is the transition probability. This probability for the rotational branching 
between two vibrational levels can be determined by the Honl-London factors. These 
factors sum to the value 2J +1 so in order to normalize the factor it must be divided by 2J 
+1. This will yield the correct branching probability. The Honl-London factors yield the 
probability of a rotational branch for a given v'-v" transition. To consider a transition on 
the specific v" level another factor Ay'y" must be included. The factor Ayy becomes 


Ajyry" = Ayty" at : (4.28) 


Equation 4.24 can now be rewritten as 


E(n'v'N'J'p'jn"v"N"J"p") = inva Se ey ; (4.29) 


1 on oe 
aAv'y 
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The single-scatter albedo @y'y" of a fluorescent band is defined as 


Oyiyihie- Eee an) 


dAyiy" 
V 


The factor Ay'y” is the transition probability of a molecular band which is related to the band 
strength by the relation shown below: 


64r4v3yiy0 
A 1 = 2 tlt 
viv the dy" Re“ qy'v". (4.31) 


In Equation 4.31 Vy'y"' is the frequency of the band, dy' is the degeneracy of the upper 
level and Re? qy'y" is a close approximation of the band strength Sy'y". These factors 


allow Equation 4.28 to be expressed as 


Svty": = ———_——._... (4.32) 


Equations 4.30 through 4.32 can now be substituted into Equation 4.29 to yield the volume 
emission rate. 

Production for forbidden transitions conforms to the required calculations 
for electron impact excitation. The relative speed and abundance of electrons in the 
ionosphere results in collisions between electrons and neutral molecules being the dominant 


production source for forbidden transitions. Excitation by electron impact becomes the 
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dominant source of the excited state for the transition. The source equation is now highly 
dependent on the electron flux, as well as the electrons cross section for excitation. The 


emission rate as compiled by Barth (1965) is given by 


Emax 


i lleoe aU chi) (4.33) 


where Ig is the particle flux in electrons cm-2 sec-!. Both, the electron particle flux and the 
electron cross section for excitation of an electronic state Og are energy dependent. The 
Franck-Condon factor for excitation qoy', single-scattering albedo @y'y", and molecular 
column density Ng, adjusted by the Chapman function ch(9), remain consistent with 
those for the photoexcitation case. It is necessary to integrate over all available electron 
energies when a range of electron energies is present. This is the usual case in the 
ionosphere. 

The methods of this section have been used to compile the synthetic spectra 
to be illustrated in Chapter V. Each spectrum is highly dependent on the unique properties 
of each molecule involved. Separate algorithms have been used for each emission feature 
in Chapter V. The actual calculation requires that a volume emission rate for each line in 
the band be calculated. The contribution from each bin is then summed to produce the 
spectrum. An example of this process is shown in Figure 4.10. In this figure the (1,0) 
gamma band of Nitric Oxide is shown for a wavelength bin size of .03 A. The individual 
lines are distinct. The spectra in Chapter V have a lower resolution comparable to 
instrument resolving capability. As a result individual rotational lines are not 


distinguishable. 
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Figure 4.10 Synthetic Spectrum for (1,0) Gamma Band Emissions 


(Cleary, 1985) 
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V. DAY GLOW EMISSION FEATURES IN THE 1800-3400 A 
WAVELENGTH RANGE 


This chapter is concerned with identifying observable emission features in the 
wavelength range of the instrument. N2, O2, O and NO all have emissions in the 1800- 
3400 A window. The transitions result from the interactions presented in Chapter IV. 
There are additional emission bands that occur in the window besides those that will be 
described. These are of such a low intensity at the altitudes of interest that they are not 
measurable. The bands and lines presented in this chapter are primarily the ones that are 
predicted to be observed, during the flight. Computer simulations of synthetic spectra have 
been created by McCoy (1981), Cleary (1985), and Conway (1982). 

The observed intensity will be limited by the detector efficiency. The instrument 
detector efficiency is limited at the long end of the experimental wavelength range by 
photocathode efficiency. Figure 5.01 shows this efficiency over the wavelength range of 
the instrument. The solid curve in the figure demonstrates that detector efficiency drops off 
rapidly past 3000A . This efficiency drop will become a significant consideration when 


attempting to identify emission features of longer than 3000A that are observable. 
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Figure 5.01 Photocathode Efficiency 


A. MOLECULAR NITROGEN 


Molecular nitrogen has three observable band systems: the Vegard-Kaplan (A2Z,+ 
- XI Xgt), 2nd positive (C3My - B31Ig), and the Lyman-Birge-Hopfield (a! Tg - X!T1g) 
bands. The Birge-Hopfield (b!Nly - X1Xt,) band was also investigated as a possible 
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contributor but its intensity in the experimental window was found to be negligible. The 
potential energy curves for N2 are shown in Figure 5.02. References to specific energy 
levels of the Nz molecule will be from this figure. In addition, Figure 5.03 is a schematic 
energy level diagram of N2. This figure illustrates the major transitions of molecular 
nitrogen. It also displays the relative energy differences for the emission bands. 

It is of note that molecular nitrogen is the dominant component of the atmosphere at 
the lower end of the experimental altitude range. However, at around 175 km atomic 
oxygen becomes the dominant species (see Figure 2.08). Since N2 is the dominant 
species, a reaction involving N2 may be observable although it has a relatively low 
emission probability. Additionally, the relative importance of N2 is strongest at lower 
altitudes and will tend to drop as altitude increases unless some other factor arises such as 
quenching at lower altitudes. The N2 emissions in the ultraviolet tend to be weak because 
they are a result of forbidden transitions subject to quenching, or a combination of both 


factors. 
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Figure 5.02 Nz Potential Energy Curve 


(Journal of Physical Chemistry Reference Data, 1977) 
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Figure 5.03 N2 Energy Level Diagram (Sharp, 1970) 
ds Vegard-Kaplan System 


The V- K system represents the A3Z,,+ - X!y ot transition, and has a range 


of 1250 to 5325A (Krupenie, 1977). In 1932, Vegard first discovered the system by 


bombarding solid nitrogen with electrons. Later in 1934, Kaplan succeeded in obtaining 
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V-K bands with gaseous nitrogen. Kaplan also adjusted the wavelengths of the lines 
proposed by Vegard. (Krupenie, 1977) 

Barth (1965) identified both fluorescence and electron processes as 
production sources for the V-K emissions. However, Sharp (1970) has determined that 
the sources for the upper state are almost exclusively due to electron impact. The upper 
state ASZ+, is populated as a result of direct excitation by electron impact, as well as 
cascading from higher states (Broadfoot and Hunten,1964). Dominant sources of cascade 
are directly from the BIT, state and indirectly from the C3I]y state. It is of note that the 
relative population of the vibrational levels depends on the excitation source. Table 5.01 
has been computed to demonstrate the relative distribution of vibrational levels as a result of 
each case. The first row considers the distribution for direct excitation to the A3Z+, state. 
The second row considers the relative distribution of vibrational states only for cascade 
from higher states. While direct excitation results in a concentration at higher vibrational 
levels, cascade results in a concentration at the lower vibrational levels. Table 5.01 is in 
non-normalized units where a larger number indicates a higher relative population of the 


entire A5Z+, state. 


TABLE 5.01 
RELATIVE POPULATIONS OF THE A32+y STATE 
(Ahmed, 1969) 
Vv: 4 l 
direct .15 1.0 2.75 6.0 10.25 15.250 19:5 2273 24°75 5 25mee2 
cascade .95 1.0 .78 .57  .35 19 09 
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This scheme is further complicated by the proposal by Campbell and Thrush (1967) and 
later by Gilmore (1969) that inverse first positive transitions act as non-radiative loss 
mechanisms (see Figure 5.03). It has been suggested by Torr and Torr (1982) that the 
complex population mechanisms for the V-K system require further data. This experiment 
may be able to provide some useful data at the upper end of the experimental altitude range. 

Because the upper state of the Vegard-Kaplan state is metastable, quenching 
rather than emission will dominate. This is particularly true at the low altitude end of the 
experiment. It has been shown that the Vegard-Kaplan band is a weak feature in the 
dayglow although it is quite bright in the aurora (Torr and Torr, 1982). The primary 
reason for this appears to be quenching. Shemansky (1969) determined that the upper state 
possessed a mean lifetime of 1.8 seconds. The atmospheric collision frequency is about 
2.25 per second at 120 km ( CIRA model 1965, in Sharp, 1970). Therefore at low 
altitudes there is a likelihood of four collisions prior to an emission. Quenching may result 
from collisions between N2 and neutral atoms or molecules. Candidates include N2 O2, 
O, NO, and N as potential quenching agents (Torr and Torr, 1982). Sharp (1970) has 
identified O as the dominant agent. The number density for all species drops by two orders 
of magnitude between 100 and 200 km (see Figure 2.08). Collision frequency can be 
expected to drop off proportionally. Thus, the Vegard-Kaplan bands will assume an 
increasingly significant role as altitude increases. 

Conway (1985) has identified the (0,9) band at 3353A as having an 
intensity of 2 kR at an altitude of 150 km. This is a detectable amount although it is well 
past the 3000A threshold of reduced detector efficiency. It is necessary to investigate other 
vibrational bands in the V-K system to determine whether or not this band system will be 
observable. There is a convenient factor to distinguish the relative intensities of emission 


bands for a particular band system. This factor is a combination of the Frank-Condon 
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factor for excitation and the albedo for scattering (Q- omega). This is interpreted as that 
portion of the total excitation to a particular electronic state that is excited to a specific v’ 
state and emitted back down to a particular v" state. The Q.omega factor for the (0,9) 
feature is 2.47x10-5. Numerous emission features below the 3000A cutoff have similar or 
even higher Q-omega factors (Barth, 1965). The highest is the (5,1) band at 1841A, with 
Q-omega equal to 1.48x10-2. This would support the contention that N2 V-K should be 
observable in measurable quantities during this experiment. 

The V-K transition is spin forbidden but remains electric dipole radiation 
(Conway, 1980). It violates the As = 0 requirement as the spin 3 to spin | states, are not 
allowed. However,the 3Z,+ to Dot transition wul occur despite it being forbidden. 
Herzberg (1950) states that the spin orbit interaction causes the !Z state to "borrow" some 
of the characteristics of the allowed 31g state, or the 3 state "borrows" some of the 
characteristics of the !II state. Four branches will result. The synthetic Vegard-Kaplan 


spectrum between 1800-3400 A is shown in Figure 5.04. 
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Figure 5.04 N2 Vegard-Kaplan Bands 


De. N2 Second Positive System 

The second positive system represents the eine BIg transition. It 
ranges from 2680 to 5460A (Krupenie, 1977). The Second Positive system was first 
observed in 1869 (Krupenie, 1977). It appears in most nitrogen sources including pure 
nitrogen and air. It was observed in the dayglow as early as 1945 (Pearse and Gaydon, 
1963). Numerous observations include a recent experiment by Conway (1980) in 1980. 

Conway and Christensen (1985) reported an observed intensity of about 
10kR for the (0,0) feature at 3370A during solar maximum conditions. It appears that this 
is the brightest feature in the experimental wavelength region. Other features in the region 


are less intense. The Q-omega factor for the (0,0) band is .280 . The highest factor below 


61 


3000A is the (2,0) with a value of .021 (Barth, 1965). This is an order of magnitude lower 
than the (0,0) band and is an indication that the 2nd-Positive band will make little 
contribution to the observed emissions in the experimental wavelength region. 

The intensity can be calculated using the following method. From Equation 
4.33, the intensity due to electron impact excitation 1s equal to the product of the total 
electron impact cross-section, slant column density, and Q-omega factors for the band. 
The calculation used the electron impact cross-sections for the C311, state which is shown 
in Figure 5.05. The cross-section at each energy can be multiplied by the electron flux at 
that energy. The electron flux has been calculated in a simulation by Strickland and Meier 
(1982) using an F10.7 equal to 206, a solar zenith angle of 56 degrees, and a temperature 
at infinity of 1OOOK. The data points on the figure indicate separate results by E. Oran 
(Oran, 1978). The electron flux at each energy can be multiplied by the cross-section at 
each energy the result is shown in Figure 5.07. The integration over the electron energy 
range shown in Figure 5.07 will yield the total electron excitation cross-section for each 
State. This cross-section multiplied by the Q-omega factor and a representative column 
density at a given altitude will give the expected observable intensity in Rayleighs. The 
slant column density is computed by adjusting the vertical column density by the chapman 
function for looking across the limb. This has been done for the values given in Table 


5.02. 
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TABLE 5.02 


INTENSITY FACTORS FOR N2(C3My - B3Tg) TRANSITION 
Emax 
Jor Ip dE = 8.25x10°9 


dov'@y'y" (Q-omega) = .28 (Barth, 1965) 
N =1.8x10!8 mol cm-2{ for Nz n(150km) = 1.0x10!9mol cm°3, ting = 1800°K, 
ch(30,90°)} 
4nI =4.2kR 


This calculated intensity is a reasonable approximation to the measured value of 10kR, 


considering that both N2 column densities and electron fluxes are approximated. 
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Electron Impact Cross—Sections vs. Electron Energy 


i, My (CT ) 


Hlectron Impact Cross—Section (cm2) 





i 10 100 1000 10000 
Electron Energy (eV) 


Figure 5.05 Electron Impact Cross-Section vs. Electron energy N2(C3Iy) 
and O2(A2Z+*q) states (Cartwright et al, 1977 and Borst, 1972). 
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Figure 5.06 Electron Energy ys Electron Flux (Strickland & Meier,1982) 
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Figure 5.07 Electron Energy vs. (Cross-section x Electron Intensity) 
Because of the the AZ = 0 selection rule discussed in Chapter V, three sub- 

bands will result for Hund's case (a) 

>Tig - 310 

311, - 31 

ahd olubs 
(Herzberg, 1950). Each sub-band has a strong P and R branch. The Q branch is a weak 
branch except for the 3Ig - 3IIg case. This results in triple band heads which are degraded 
to the blue. The Physical Chemistry Reference Data Manual (Krupenie, 1977) has 
identified fifty of the bands. When J is small, the transition conforms to Hund's case (a) 
(see Figure 4.07). As J becomes larger, the coupling approaches that of Hund's case (b) 
(see Figure 4.08). This transition exhibits some A-type doubling doubling, as required by 
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Hund coupling for a 37 - 31 transition, The second positive synthetic spectrum that will 


contribute to the experimental window is shown in Figure 5.08. 
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Figure 5.08 Nz Second Positive Bands 
3% Lyman-Birge-Hopfield Band 
The Lyman-Birge-Hopfield bands (a! TI], - X!Z*,) have a wavelength 
range from 1160 to 3020 A (Barth, 1965). The emission bands were first observed in 
1928 by Birge and Hopfield (Pearse and Gaydon, 1963). They are readily visible as both 
emission and absorption features in the vacum ultraviolet. 
The transition is forbidden as electric dipole radiation Emissions will result 


from electric quadrupole and magnetic dipole transitions. Furthermore, the upper state is 
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forbidden from the ground state (Herzberg, 1950). The excitation will occur primarily as a 
result of electron bombardment. 

LBH emissions have been used as an accurate measure of the photoelectron 
flux (Conway, 1982 ). Conway observed dayglow intensities approaching 800 R at 160 
km for the (2,0) feature at 1383 A. The Q-omega factor for the (2,0) band is 4.63x10-2. 
The (2,8) band at 1837A has a Q-omega factor of 1.87x10-2. The (3,9) at 1854A, (3,10) 
at 1928A, (4,10) at 1870 and (4,11) at 1945A emission bands all have similar Q-omega 
factors; on the order of 10-2 [see Barth (1965) for the complete table electron impact cross- 
sections]. This would indicate that LBH emissions should be visible in measurable 
intensities in the wavelength region. These features are at the short wavelength end of the 
detector which has a relatively low efficiency . Nevertheless, it would seem that the 
Lyman-Birge-Hopfield Bands should contribute to the observed experimental intensity. 
The peak observed intensity for all LBH emission features should be in the range of 150 to 
200 km in altitude if atmospheric conditions during Conway's observations remain 
generally in effect. 

The a! TI, - X!£*, transition violates the g ¢ u selection mule given in Table 
5.1. The transition is allowed by the special selection rules of Table 5.2. The AJ selection 


rule for quadrupole radiation becomes 


AJ =0, +1, + 2 (Herzberg, 1950). 


The resultant branch structure is some what complicated. It is best presented graphically as 


is shown in Figure 5.09. 
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Figure 5.09 II, - 1£+, Branches for Magnetic Dipole and Quadrupole 
Bands (Herzberg, 1950) 
The branch structure of the LBH system will produce the lines shown in Figure 5.10 in the 


experimental window. 
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Figure 5.10 N2 LBH Bands (Cleary, 1985) 


69 


B. MOLECULAR OXYGEN 

There are two observable bands of O2 that may contribute to the experimental 
window. These are the Herzberg I, and Schumann-Runge (SR) bands. The first of these 
is forbidden, while the second is an allowed transition. The potential energy level diagram 
for O2 is shown in Figure 5.11. The schematic energy level diagram for O2 is shown in 
Figure 5.12. This figure demonstrates the major transitions of molecular Oxygen. It also 
displays the relative energy differences for the emission bands. The Schumann-Runge 
bands have an easily accessible threshold energy for photo dissociation which minimizes 
the observed intensity. The energy curve for the B3Z-y is of note. This is the upper state 
for the Schumann-Runge transition. O2 in this state will dissociate into O(3P) + OUD) 
atoms. It may also predisscociate to the O(3P) + O(3P) atoms. The upper state for the 


Herzberg system (A3Z+y) may also result in the O(7P) + O(3P) atoms. 
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Figure 5.11 O2 Potential Energy Diagram (from Gilmore, 1965) 
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bes Schumann-Runge Bands 


The Schumann-Runge system (B?X-y - X3Z-y) ranges from 1750 to 5350 
A. With over 100 single-headed red degraded bands, the Schumann-Runge (SR) is a very 
extensively studied Oxygen system. It was identified by Schumann as early as 1903 
(JPCRD, 1972 ).The large number of lines arises from the greater internuclear distance of 
the upper electronic state over the lower state. This results in a large number of accessible 


lower vibrational levels for each upper vibrational level. 


The B32-,, state has a strong absorption cross section However as 





mentioned above, this state lies very close to the Oxygen dissociation continuum. The 
ionization and dissociation continuum are shown with respect to the absorption cross 


section in Figure 5.13. 
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Figure 5.13 Absorption Cross-Section for O2, Showing 
Continua, and SR Bands (Banks, and Kockarts, 1973) 

For v levels greater than three, the O2 molecule will predissociate rather 
than fluoresce (Banks and Kockart, 1973 ). This means only the v'=0,1,2,3 states will 
contribute to the emission spectrum. In fact the Schumann-Runge system is actually the 
major source of atomic oxygen. 

An estimate of the SR intensity can be made by comparing its g factor and 
the O» slant column density with the g factor and slant column density of known 
emissions. Barth (1965) lists the number of photons scattered per second per atom (g 
factors) for numerous atmospheric constituents. These g factors multiplied by the slant 
column density determine the relative strength of an allowed transition. The g factem for 
the first four v' states for the Schumann-Runge bands are extremely small. The largest 
value in the experimental wavelength range is 3.59x10-!! for the (3,7) band at 2442 A. To 
put this in perspective a bright allowed feature such as the Nitric Oxide (1,0) gamma band, 


has a g factor of 7.7x10-®. This implies that an NO Gamma transition is much more likely 
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than an O7 Schumann-Runge transition. However, a contrasting value from Figure 2.08 
shows that the number density for O2 is about five orders of magnitude higher than NO. 
This would seem to imply that the strongest O2 SR bands should be of comparable 
intensity to NO gamma bands. The values used to calculate the (3,7) SR band at an altitude 


of 150km for active solar conditions 1s provided in Table 5.03. 


TABLE 5.03 
O2_SR3,.7BAND INTENSITY CALCULATIONS 
g = 8.25x10°9 
N =1.8x10!8 mol cm-2{ for O2 n(150km) = 1.0x109 mol cm°3, ting = 1800°K, 
ch(30,90°)} 
4nI =4.2kR 


The Oxygen Schumann-Runge band system was not identified by Sharp (1986) in the 
thermosphere between 2000 and 8000 A. No other experimental observation of this band 
has been found. However, it appears that emissions do exist in the experimental 
wavelength region, and they may contribute to the observed intensity. 

The 3Z-, state is an allowed state for photon excitation from the X3E-5 
ground state. Once in the excited state the 3Z"y - to 3Z*, transition produces both P andR 
branches. Each branch has three components for a total of six main branches for each 
band. Additionally there are six satellite branches where AJ # AK. Splitting will become 
significant only when K becomes large (Herzberg, 1950 ). The O2 molecule has zero 
nuclear spin. This condition requires the molecule to follow Bose statistics (Herzberg, 
1950 ). This requirement produces a band structure where every alternate triplet in the P 


and R branches appears. The band system is red degraded. 


ié 


The synthetic spectrum for the SR bands that will fall in the window is 
shown in Figure 5.14. The synthetic spectrum shown is for all the v'=0,1,2,3. For v'>3 


the molecule will predissociate rather than fluoresce. 
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Figure 5.14 O2 Schumann-Runge Bands 
Ze Herzberg Bands 
The Herzberg transition (A%Zty - X3X°g) appears to range from 2590 to 
4880 A (JPCRD, 1972 ). The band system appears to have been discovered rather recently 
by Herzberg, in 1952 (Pearse and Gaydon, 1963). The band system 1s a prominent feature 


in the nightglow, where the excitation source is inverse predissociation (Barth, 1965 ). 
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As a forbidden transition, Herzberg bands will result primarily from 
photoelectron impact. Its dayglow intensity will be less than a corresponding transition in 
the Schumann-Runge Band. The 2+ — Zt or 2 <2 - selection rule is the selection rule 
for this transition. 

The upper state is subject to dissociation for incident energies greater than 
5.1eV. It corresponds to a wavelength of 2420A. This means that at any altitude, 
emissions as a result of the A3Z+y - X3Z-,g transition will be longward of 2420 A. The 
Herzberg dissociation peaks at 40 km and is closely linked with ozone production 


(Chamberlain, 1987). The A3Z+, state dissociation reaction is shown below: 
O2+e"(E> 5.leV) => OP) + OP) +e. 


The calculation to determine the intensity of the Herzberg bands follows 
from that developed for the N2 Second Positive bands. Figures 5.05, 5.06, and 5.07 are 
used to determine the intensity. The calculations were done for active solar conditions, 
looking across the limb. The brightest band identified is the (12,1) band at 2528A. The 
factors used to determine the observed intensity are shown in Table 5.04. 


TABLE 5.04 
INTENSITY FACTORS FOR O2(Herzberg) TRANSITION 
Emax 


doz Ip dE = 1.21x10-7 


dov'@y'y" (Q-omega) = 1.98x10-3 (Barth, 1965) 
N =1.8x10!7 mol cm-2{ for O2 n(150km) = 1.0x109 mol cm°3, ting = 1800°K, 
ch(30,90°) } 
4nlI =45R 
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The calculated intensity of 45R is a marginal value as an observable value in the wavelength 
region. The Herzberg bands will contribute to the total observed intensity, however 
separate emission features may not be distinguishable. The O2 Herzberg spectrum that 


contributes to the experimental window is shown in Figure 5.15. 
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Figure 5.15 O2 Herzberg Bands 


OF NITRIC OXIDE 

Nitric Oxide has three band systems in the experimental window. The Gamma 
band (A2Z+ - X2I1), Delta band (C2II - X2I1), and Epsilon bands (D2Z+- XII). The 
potential energy diagram for NO is shown in Figure 5.16. A schematic energy level 


diagram, which illustrates the major NO transitions, is shown in Figure 5.17. 
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Although NO has a significantly lower density compared with N2, O2, and O, the NO 
bands are very intense. As shown in Figure 2.08, at the lower end of the experimental 
altitude it is approximately five orders of magnitude less than the N2 and O> densities. 
However, all of the the NO band systems are allowed transitions. The importance of nitric 
oxide as an emission source of the airglow becomes as significant as a corresponding 
forbidden transition of a more abundant species. The importance of NO diminishes with 
altitude beyond its peak number density at about 105km. A fourth system, the 

Beta band (B2I - X21), also has transitions in the experimental window. However, the 
spatial overlap for this system results in intensities four to five orders of magnitude less 


than the other bands. It is not a significant contributor to the window. 
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Figure 5.16 NO Potential Energy Diagram (Cleary, 1985) 
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Figure 5.17 NO Energy Level Diagram 
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lies Gamma Band 

The Gamma Band (A2=+ - X2IT) is the most intense of the NO transitions 
in the experimental wavelength range. The bands range from 1873 to 6126 A (Barth, 
1965). The band system was observed as early as 1917 by Strutt (Pearse and Gaydon, 
1963). It commonly occurs in flames containing nitric oxide as well as shock-treated air. 
The (1,0) Gamma band, at 2149 A, is the most intense feature of Nitric Oxide in the 
ionosphere. The Gamma bands were not observed in the ionosphere until 1964 by Barth 
(1965). 

The primary excitation mechanism for the upper state is photoexcitation. 


Other sources to the A2Z+ state include cascade from the C2Z+ state(Sharp, 1986) and the 


chemical reaction shown below: 


NO(B2I; v=7) + M = NO(A2Z+; v=3) + M(Young & Sharpless, 1962) 


These sources result in emission, where the NO molecule returns to the XI ground state . 
The upper state has a lifetime of 2x10-7 seconds (Callear and Pillings, 1970) which insures 
emissions rather than collisions will dominate activities in the state. 

Measurements of the most intense Gamma bands are limited by self- 
absorption. Self-absorption has been observed by Cleary (1985) below 125 km for the 
(1,0) band at 2149A, and below 115 km for the (2,0) band at 2047A. These two bands are 
the strongest emission features of the gamma band system. Self-absorption is proportional 
to the column density times the oscillator strength for a band. This is why self-absorption 
is most significant for the strongest bands. 

Nitric Oxide Gamma bands have been observed recently by Cleary (1986) 


in 1982, during active solar conditions. His observations yielded observed intensities of 
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about 30kR for the (1,0) band at an altitude of 105 km. This intensity had diminished to 
2kR by an altitude of 190 km. Another measurement of the NO (1,0) gamma band at 2150 
A by Thomas (1978) was conducted during solar minimum conditions. His measurements 
indicate a measured slant intensity of 10 kR at 100 km, diminishing to less than 100R at 
300 km. These figures support the prediction that the NO Gamma band will be a 
significant contributor to the experimental window, especially at low altitudes. 

The upper state A2Z+, conforms to Hund's case (b) while the ground state 
is an intermediate case of (a) and (b). The 2 - 2I1 transition results in doublet splitting for 
both the upper and lower states. This produces a total of four band heads. Additionally, 
the AJ =0, + 1 selection rule will result in P,Q, and R branches. Therefore the entire band 
will contain twelve branches. Bands are degraded to the blue (Herzberg, 1950). The 


emission spectrum of the NO Gamma band emission spectra is shown in Figure 5.18. 
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Figure 5.18 NO Gamma Bands (Cleary, 1985) 
Ds Epsilon Bands 
The NO Epsilon bands (D2Z*- X2I1) are the same type of electronic 

transition, Z+- 2II, as the gamma bands. The band system ranges appears to be most 
intense in the ultra-violet although precise wavelength ranges have not been identified. 
This emission system is one of the lesser known in the NO family, although it was 
observed by Guillery in 1927 (Pearse and Gaydon, 1963). Separate emission features in 
the ionosphere were identified by Cleary (1986). 

The emission is a similar resonance scattering as the gamma bands. This means 


that an excitation from the ground state to the excited state leads to a relaxation back to the 
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ground state. The exclusive excitation mechanism for the D2Z*+ state is photoexcitation. 
(McCoy, 1981 ). 

Cleary (1986) measured a peak intensity approaching 700 R/A for the (0,1) band at 
1950 A at an altitude of 105 km. It appears that this band as well as the slightly weaker 
(0,2) at 2020A and (0,3) at 2090A features should be observable in the experiment. 

Once again each band has four band heads, which yields the same twelve branch 
structure as the gamma bands. As shown in the potential energy diagram, an epsilon 
transition will be generally of higher energy than a corresponding gamma transition. The 
NO epsilon bands that will contribute to the experimental wavelength region are shown in 


Figure 5.19. 
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Figure 5.19 NO Epsilon Bands (Cleary, 1985) 
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Jy NO Delta Bands 
The Delta Bands (C2II - X2I1) are also an allowed transition for the 
211 - 211 case. Delta bands range from 1621 to 4288 A. Early observations of this band 

system were made by Knauss in 1928, and Schmid in 1930 (Pearse and Gaydon, 1963). 

The Delta excitation mechanism is photoexcitation similar to the Epsilon and 
Gamma bands. Delta bands do not demonstrate a fluorescence efficiency of unity, 
however. The accepted efficiency factor is .25 (Cleary 1986 ). This low efficiency is 
required because of predissociation. 

Cleary (1986) also observed delta bands during his 1982 experiment . 
Sharp (1986) has not identified any more recent observations. The strongest feature 
observed by Cleary was the (0,1) band at 1980 A. He measured a peak intensity of about 
150 R/A which will be difficult to distinguish in the experimental wavelength region. 

Once again each band has four band heads, which yields the same twelve 
branch structure as the gamma bands. As shown in the potential energy diagram, a delta 
transition will be generally of higher energy than a corresponding gamma transition, 
although it is of a slightly lower energy than the epsilon bands. The NO delta bands that 


will contribute to the experimental wavelength region are shown in Figure 5.20. 
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Figure 5.20 NO Delta Bands (Cleary, 1985) 
Oe ATOMIC OXYGEN 
There is one atomic oxygen line that occurs in the experimental window. This is 
the 2972 A line. The line is a result of the (O(!S) - O(P)) transition. It is a well known 


forbidden emission of the oxygen atom. 


This line is represented by the transition shown below. 


O(S) = O(3P) + hv(2972A) 
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This transition is forbidden as it violates AS = 0 requirement for atomic transitions. 
Although forbidden, this transition has been reported by Sharp at an average intensity of 
80R. This is a marginal observed intensity. Furthermore, this line is near the 3000A 
instrument efficiency cutoff which will tend to further limit the observed intensity. 
However, the relative density rises steadily throughout the experimental altitude range. The 
2972A line will be swamped by stronger emissions at low altitudes in the experimental 
region but could be distinguishable at higher altitudes. This will be especially true when O 
becomes the dominant constituent above 175 km. 

The upper state for this transition, the O('S) state, is the source of two branches. 


These two branches are shown in the energy level diagram of Figure 5.21. 
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Figure 5.21 O Energy Level Diagram 
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It is clear from Figure 5.21 that an O atom in the |S state can either emit a 2972 or 5577 A 
photon. Furthermore, since the branching ratio between the 5577 and 2972 A lines has 
been determined a measurement of the 2972A feature can be used to determine source terms 
for the O('S) state. The sources for the O('S) state are shown in Table 5.05. 


TABLE 5.05 


Oot+N 
N2(A3Zyt) + O 
These sources include N? in the upper state for Vegard-Kaplan bands. This relationship 
could be used to assist in determining the quenching rate for the A7Z+, state of Nz. Cross- 
sections and efficiencies for these source terms are given by Sharps (1986). 
The 2972 A line should be observable in the experimental window. It is near the 
limit of optimum detector efficiency, but below the 3000 A threshold of reduced efficiency. 


It appears that this line may also be useful in determining N2 molecular densities as well. 
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VI. CONCLUSION 


Pay: SUMMARY OF FINDINGS 

Identification of emissions in the wavelength range of the MUSTANG experiment 
was accomplished using experimental results to support calculated emission feature 
intensity. All emission bands that will contribute to the experimental spectrum, as well as 
those features that will be distinguishable, have been tentatively identified. These results 
confirm the contention that the instrument wavelength range allows observation of all 
midultraviolet airglow emission bands. The extended wavelength and altitude range for 
this experiment allows a unique experimental opportunity to observe all major neutral 
atmospheric constituents in one spectrum. 

The major midultraviolet emission bands are shown in Table 6.01. This table also 
identifies the brightest emission feature, its wavelength, and whether or not the band will 
be distinguishable. The last column provides special circumstances for observations of that 
band system. This column considers the most intense feature and accounts for the 3000A 
instrument efficiency cutoff, where appropriate. It can be seen in Table 6.01 that the 
spectrum is quite complex. Analysis of experimental results will require unfolding each 


constituent spectrum from the combined observed spectra. 
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TABLE 6.01 


MUSTANG EXPERIMENT: OBSERVABLE EMISSIONS 


Species System Band ACA) Observable 

N2 VK (5,1) 1841 above 200km 

2P (2,0) 2973 marginal due low Q-@ 
(1,0) 3157 marginal due to inst. eff 

LBH (3,9) 1854 yes 

O2 SR (3,7) 2442 yes 
H-I (12,1) 2528 marginal 

NO Y (1,0) 2149 yes 
E (0,1) 1950 yes 
) (0,1) 1980 marginal 

O line N/A 2972 above 175km 


B. TOPICS FOR FURTHER INVESTIGATION 

This preliminary identification of the predicted experimental spectrum for the 
MUSTANG experiment must initially be confirmed by data analysis to identify those 
emission features in Table 6.01. This data will then be analyzed to identify N2 density in 
the E, Fl, and F2 regions along with NO density tn the E region. Further analysis along 
with data from the HIRAAS instrument will be used to identify the solar EUV flux and 
photoelectron flux. These results will contribute significantly towards meeting the DOD 
goal of identifying the ionospheric electron density. The experimental wavelength range 
allows observation of all major neutral contributors and a self consistent model of the 


electron density may be possible when analysis of all experimental results is complete. 
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